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ABSTRACT 

We present a theoretical study of the deuterated species detectability in various types of ex- 
tragalactic star-forming regions based on our predictions of chemical abundances. This work is 
motivated by the past and current attempts at observing deuterated species in external galaxies 
such as NGC 253, IC 342 and the LMC. Here, we investigate the influence of the density, the 
temperature, the FUV radiation field, the cosmic ray ionisation, and the metallicity on the frac- 
tional abundances and D/H abundance ratios of about 20 deuterated species. Without modelling 
any particular source, we determined how the deuterium chemistry behaves in different physical 
environments such as starburst, cosmic-rays enhanced environments, low metallicity and high 
redshift galaxies. In general, our predicted column densities seem in good agreement with those 
derived from the current limited dataset of observations in external galaxies. We provide, for the 
first time, a list of key deuterated species whose abundances are high enough to be possibly de- 
tectable by the Atacama Large Millimeter Array (ALMA) and Herschel, as a function of galactic 
nuclear activity and redshift. 

Subject headings: Galaxies: nuclei-ISM - Submillimeter - ISM: molecules - Astro-chemistry - Methods: 
numerical 



1. Introduction 

The past few years have seen a boom in stud- 
ies of molecular deuteration in galactic low-mass 
and high-mass star forming regions, triggered 
by the discovery of a large fraction of mono-, 
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chemistry, especially in highly embedded regions, 
is important because in such environments, CO 
is depleted and deuterated species are thus driv- 
ing the chemistry. In embedded regions in our 
Galaxy, a large deuterium fractionation is ob- 



served, well above t he elemental a bundance ratio 
D/H of 1.5 x 10~ 5 (jOliveiralbOOlh High molecu- 
lar deuteration can occur by two different mech- 
anisms. One involves gas phase chemistry and 
ion-molecule deuterium exch ange reactions taking 
place at low temperatures ( Watsonl Il980l ). The 
second mecha nism is based on grain chemistry 
(|Tielenslll983h . 

In this work, we investigate deuteration in ex- 
tragalactic star-forming environments, where the 
physical conditions can be drastically different 
from what we see in our own Galaxy. To date, very 
few detections of deuterated species have been 
attempted in extragala c tic environments (e.g. 
Mauersberger et ah 1995 : Martin et al. 2006). In 
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a previous paper, Bavet et al.l (|2008l ) investigated 
the non-deuterated chemistry in a variety of ex- 
tragalactic environments and obtained theoret- 
ical predictions for the warmer and dense gas, 
traced by non-deuterated species. Here we im- 
plement the mono-deuterated species (hereafter 
ca lled D-specie s ) netw ork to our existing one used 
in iBavet et al.l (|2008l ). in order to investigate the 
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differences in the D- fractionation and the D/H 
ratios in regions of star-formation (i.e. gas with 
n(H2)> 10 5 cm -3 ) in various types of galaxies (i.e 
starburst, cosmic-rays enhanced environments, 
low metallicity and high redshift). We also aim 
at guiding observers with potentially detectable 
deuterated tracers of so called dense star-forming 
molecular gas in external galaxies when using Hcr- 
schcl or the forthcoming ALMA, as based on the 
strength of their chemical abundances. 

The paper is divided in five sections; Section [2] 
describes the model used and presents the choice of 
model parameters used for mimicking various ex- 
tragalactic environments. Section [3] presents the 
D-species fractional abundances and D/H abun- 
dance ratios obtained in the various gas compo- 
nents and environments investigated. Finally, in 
Section HI we discuss the results and conclude. 

2. Model descriptions 

We have used the same chemica l model as 
the one used in iBavet et al. (2008). but with 
an extended chemical network to include mono- 
deuterated sp ecies. The chemical cod e, U CL_Chem, 
i s pres ented in Viti fc Williams ( 19991) and Viti et al 
(|2004l ). We refer the reader to these papers for a 
full description of the code. 

To summarise, UCL_Chem model is a time 
and depth depen dent chemical model. As in 
Bavet et"all ((2008|), we first model the collapse of 
a 10K core (Phase I); we then follow the chemi- 
cal evolution of the region once the stars are born 
(Phase II). The temperature in both Phases is the 
typical average for starless and hot cores, respec- 
tively, since here we are considering an ensemble 



of cores within a large region of space. I Bavet et al 



( 2010bl ) do in fact explore in great detail the influ- 
ence of changes in temperature and density within 
a single hot core when subject to various environ- 
mental conditions. They find that thermal varia- 
tions due to temperature structure would flatten 
out when averaged over a large number of cores. 

The presence of an infrared source in the centre 
or in the vicinity of the region is simulated by sub- 
jecting the gas and the dust to an increase in tem- 
perature. In our models, we have assumed equal 
dust and gas temperatures. This is justified since 
in Phase II the opacity is very high (see values 
listed in Table 2]) and thus the thermal coupling 



between gas and dust is valid. 

In both phases, the chemical network, now im- 
proved with D-species, is based on more than 4000 
chemical r eactions adapted from the UMIST 2006 
database (iMillar et al.lll997t iLe Teuff et aljlioool 
IWoodall et al.l l2007h involving about 300 species 
(92 of them being deuterated) of which 62 are sur- 
face mantle species. The relevant surface reactions 
included in this work involve mainly simple hydro- 
genation, apart from the formation of the CH3CN 
which is believed to be formed via reactions in - 
volving HCN on grains (jGarrod fc Herbst! 12006). 
We note that our approach for surface reactions is 
rather simplified with r espe ct to that adopted b y 
Garrod fc HerbstJ (|2006l ) and lCazaux et alj(|2010h . 



Nevertheless, a qualitative comparison of the man- 
tle compo sition at the end of P hase I with the re- 
sults from lCazaux et al" (|2010l ) mantle show good 
agreement (c.f. H 2 and HDO abundances on 
grains as a function of the radiation field for ex- 
ample). As a first approximation and to keep the 
chemical network to a reasonable size, we include 
only the mono-D species. D2CO and its corre- 
sponding ion are the only doubly-D-species in our 
network. 

The deuterium chemical netw ork i s similar to 
that o f lRoberts fc Millar! (|2000alibh and lRoberts et al 
(|2004l ) . We assume that D atoms react in the same 
way as H atoms when reacting with other atoms 
and molecules. 

One of the outputs of the UCL_Chem is the 
fractional abundance (with respect to the total 
number of hydrogen nuclei and as a function of 
time) of gas and surface species from which we 
calculate the D/H ratios. 

The above assumptions (i.e. no doubly- 
D-species except D2CO formed on grain nor 
triply-D-species included in the code) do af- 
fect the reliability of the predictions for a cou- 
ple of species namely D 2CO, H2D" 1 " and N2D + 
(see Roberts et al.l [20031) . On the other hand, 
only mono-D-species have been observed so far 
in extragalactic environments, and the relative 
doubly-D/mono-D or triply-D/mono-D species 
abundance ratios seen in galactic sources are 
small (see e.g. CD 3 OH/CH 2 DOH = 0. 46 and 



CHD 2 OH/CH 2 DOH = 0.2 fromlParise et al.l 2004 . 
ND 2 H/NH 2 D = 0.05 from lRoueff et alj|200a and, 
D 2 CO /HDCO=0.02-0.33 from iRoueff fc Gerinl 
20031 ). 
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To investigate the deuterium chemistry in high- 
mass star-forming regions in various extragalactic 
environments, we have run a grid of 15 chemical 
models (see Tables [TJ [5] and [3]). For each type 
of galaxy (see below), we have studied the deu- 
terium chemistry coming from a moderately dense 
(n(H 2 ) — 10 5 cm~ 3 - index 1 used in tables and fig- 
ures), a dense (niH^) = 10 6 cm -3 - index 2) and a 
very dense gas (n.(i?2) = 10 7 cm -3 - index 3) com- 
ponent. This has allowed us to potentially dis- 
tinguish which gas component produces the high- 
est fractional abundances of D-molecules, in each 
type of environment. This also allowed us to de- 
termine which gas component is likely to be the 
best match to the observed extragalactic D/H ra- 
tios (see Sectj3|). 

For reproducing the various types of galaxies, 
we have assumed the following : 

The Normal Spiral (NS) case assumes 
standard parameter values (symbol "ST") 
similarly as previously performed by Bayct 
et al. (2008; see their Tables ffl and |2J) . 

For the Starburst (SB) case, we have used 
a FUV radiation field of 1000 times the stan- 
dard value, and a temperature of 500 K, in 
Phase II, instead of 300 K. These models aim 
at mimicking environments such as NGC 253 
or M 82. 

The Cosmic-ray enhanced environ- 
ments (SB+) case is represented in our 
models by a cosmic ray ionisation rate of 
100 times the standard value. This type 
of models aims at reproducing supernovae 
and massive star formation as seen in galax- 
ies such as Ar p 220 (supernovae rate of 4 
year -1 , see e.g. IParra et aDl2007l ). 

The Low metallicity (Low-met) case is 
represented, in our models, by a decrease in 
the metallicity to a fifth of the solar value 
coupled with an increase of the gas-to-dust 
mass ratio by a factor of 5. Indeed, in the 
code, quantities such as metallicity, optical 
depth (A v ) and the H2 formation rate coeffi- 
cient are coupled (see values presented in Ta- 
ble[4]for each model). In this galaxy type, we 
have also assumed an increase of FUV radi- 
ation field of 1000 times the standard value, 



and a temperature of 500 K since we are par- 
ticularly interested in sources such as IC 10, 
which also host a strong star-formation ac- 
tivity. 

Finally, the High redshift (High-z) case 
is represented by models with a metallicity 
(and its coupled quan tities) of a fifth sola r 
as measured in IC10 (jZaritskv et al.l[l994T ). 
which is a local dwarf galaxy often used 
as a good template for arche typical higher 



z galaxy populations ( e.g. iMadden et al 



1997 t iLerov et al.ll2006t lYin et alJl20l6Th 



increase of FUV radiation field by 1000 times 
the standard value; a temperature of 500 
K and an increase of the cosmic ray ioni- 
sation rate by a factor of 100 times the stan- 
dard value. Indeed, in distant objects (at 
rcdshifts greater than 0.1), it is expected 
to often have a combination of several nu- 
clear activities such as AGN, SB and s u- 
pcrnovae (see for instance Seymour 20091 or 
the studies on the q uasar APM08+27 9 lo- 
cated at z^ 4 from WeiB et al. 2007 and 



Riechers et al. 20091 ). We note that we may 
indeed underestimate the metallicity at high 
redshift in our study since opti cal observa- 
tions of quasars at z~ 6 (e.g. I Jiang et al 



20071: lJuarez et all l2009h recently revealed 
that a solar metallicity value is more likely 
than a sub-solar one for reproducing their 
observations. However, we emphasize that 
we have not tried in our present models to 
match the physical conditions for any high-z 
source in particular but that we model here 
what one can consider as an archetypical 
high redshift environment. Since we know 
that the early Universe is more populated by 
dwarf galaxies than any ot her galaxy type 



(e.g. IWhite fc Frenkl 119911 ). we adopt the 



values for IC10 as our high redshift standard. 

All the details of the parameters used in the 15 
models are summarized in Tables [T] to 0] 

3. Sensitivity of deuterated chemical abun- 
dances and D/H abundance ratios to 
variations in the physical and chemical 
parameters. 

The main aim of this work is to study the sen- 
sitivity of deuterated species in models of star- 
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forming regions to variations in physical and chem- 
ical parameters that may be characteristics of vari- 
ous galaxy types. A collection of likely detectable 
tracers as seen by the strength of their chemical 
abundances will be given in Subsect. 13.31 but first 
we simply analyse the trends the deuterium chem- 
istry follows as the environment (i.e. galaxy type) 
and the gas density vary. These trends are sum- 
marized in Figs. [T]to[6l 

3.1. Deuterated fractional abundances 
sensitivity 

3.1.1. Influence of the gas density 

When the gas density increases from n(H2) = 
10 5 cm~ 3 to n(H2) = 10 7 cm -3 , regardless of the 
galaxy type, we can divide deuterated species into 
three categories: those which are insensitive to 
density changes (e.g. HDO, NH 2 D and HDCS), 
those that increase with density (e.g. HDS, C2D 
and CH3OD) and those which show a decrease 
in their fractional abundances (e.g. DCN, DNC 
and CH2DOH). The fractional abundances of the 
D-species belonging to the first category actually 
tend to converge after ~ 10 5 yrs, making these 
species good D-tracers of dense gas (n(H 2 ) > 
10 5 cm -3 ) whatever the environments since they 
show abundances above 10 10 . On the contrary, we 
show in Figs. [T] to [3] that especially DCN and 
DNC decrease with an increase of density (e.g. 
compare Models NSi and NS3). This is no sur- 
prise since the same behavior is seen f or both HCN 



and HNC when the density increases (jAwad et al 
201(1 ). This phenomenon reinforces the conclu- 



sions of lBavet et all (j2009t ) that HCN, HNC and 
their deuterated counterparts are not particularly 
good tracers of the very dense gas in galaxies. 
HDS, C 2 D and CH 3 OD increase with density. The 
enhancements in their fractional abundance, by 
more than a factor of 5 (Figs. Q] to [3]) are partic- 
ularly interesting since observing their molecular 
emissions may allow us to discriminate between 
the different gas components in a galaxy. 

3.1.2. Influence of the FUV radiation field and 
the temperature (SB case) 

Variation in FUV or temperature do not lead to 
significant changes in abundances for species such 
as HDO, NH 2 D, DCN, D 2 CO, HDCS, CH 2 DOH 
and CH3OD (i.e. variations in chemical abun- 



dances less than a factor of 2-3). However, HDCO 
and HDS abundances decrease between a factor 
of 5 and 15 when the FUV radiation field in- 
creases. This is due to dissociations by FUV pho- 
tons. At early times (< 10 5 yrs), and at low den- 
sity (n(H 2 ) = 10 5 cm~ 3 ) the abundance of C 2 D 
increases by three orders of magnitude with the 
increase of both FUV radiation field and temper- 
ature. Hence, this species is likely to be a good 
D-tracer of excitation in galaxies since very sen- 
sitive to changes in both temperature and FUV 
radiation field. At later times, the increase of the 
C 2 D abundance is less pronounced, but still sig- 
nificant (i.e. about two orders of magnitude). On 
the other hand, DNC shows the largest decrease 
i.e. by more than one order of magnitude, when 
the FUV radiation field and the temperature are 
increased (linked as seen in Table [2|). This may 
be due to cosmic-ray-induced photodissociation, 
which has a significant temperature dependence. 

3.1.3. Influence of the cosmic ray ionisation rate 
(SB+ case) 

Increasing the cosmic ray ionisation rate (£) 
by a factor of 100 (Models SBj+ to SB 3 + - 
See Table [3] and Fig. [2]) leads to most abun- 
dances been significantly reduced by several or- 
ders of magnitude (e.g. HDCS, NH 2 D, D 2 CO, 
HD CO, CH 2 DOH HDO , DN C). As already found 
bv iBavet et al.l (120081). by Meiierink fc Spaansl 
(|2005l ). and iMeiierink etahl j2006|), the effects of 
increasing the cosmic ray ionisation rate are com- 
plex, and they generally lead to a chemistry that 
approaches steady-state more quickly (see FigEJ) . 
The general decrease of abundances seen when 
£ increases is probably due to a release of reac- 
tive ions such as ionized carbon, reacting quickly 
with many oxygen-bearing species such as HDCO, 
D 2 CO, etc. However, some species such as DCO + , 
N 2 D+ (both at low density only), DC 3 N and DCN 
(whatever the gas component) show abundances 
increased by several orders of magnitude with £. 
In the case of nitrogen-bearing species, we suspect 
that the main reservoir of nitrogen, N 2 , becomes 
ionized and gives rise to N + , a reactive species 
that promotes the formation of nitrogen-bearing 
species. This behaviour is also found for th e 
hydrogen counterparts (e.g. Bayet et al. 2008h . 
What is interesting here is that D-species seem 
to be affected in the same way as their hydro- 
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gen counterpart. Species such as N2D 4 " and DCN 
are thus likely to be good D-tracers of cosmic ray 
ionisation-excited environments since they show 
high fractional abundances and are sensitive to 
their change. In principle, we thus can use deuter- 
ation to constrain ( although the effects of high 
cosmic rays ionisation rates on the chemistry are 
very complex (Bayet et al. 2010, in prep) and 
highly dependent on initi al conditions as we ll as 
the available gas coolants ( Bavet et al. 2010al) . 



3.1.4. 



Influence of the metallicity (Low-met 
case) 



Low metallicity environments such as those 
described in Sect. [2] (see Models Low-meti to 
Low-met 3 in Table [3J, lead to few surprises for 
most species in that their abundance is sim- 
ply reduced accordingly (see Fig. [3] for HDCO, 
D 2 CO, DNC, NH 2 D and HDS). However, some 
species, such as C2D or DCN, show an increase 
in their abundances with a decrease of the metal- 
licity whereas some other species such as HDCS, 
CH2DOH, CH3OD do not seem sensitive to metal- 
licity changes. Therefore these latest species (po- 
tentially) are good D-tracers of dense gas (since 
they have high abundances) whatever the metal- 
licity (especially interesting at high redshift), 
similarly to their hydroge n counterpa r t (suc h as 
CH3OH) already n oted bv lBavet etall (|2008h and 
iRollig et all (l2007h . 



3.1.5. Influence of the above parameters coupled 
together (High-z case) 

In extreme conditions such as those used for 
mimicking high redshift environments (see Models 
High-zi to High-z 3 in Table |3] and Fig. [3j, most 
of the chemical abundances of D-species decrease: 
CH3OD, HDS disappear completely from Fig. [3] 
as compared to the normal spiral case (see Fig. [TJ; 
NH 2 D, D 2 CO, DNC, and CH 2 DOH have chemical 
abundances above the limit of detectability only 
at early time (i.e. before 1.5-3xl0 4 yrs) then they 
drop severely at later time; HDCS, HDO, HDCO, 
even if still above the limit of detectability at all 
times, decrease by more than two orders of magni- 
tude in high redshift environments as compared to 
the normal spiral case. These species nonetheless 
remain the best D-tracers of dense gas components 
in high redshift environments since very sensitive 
to the parameters changes. Only DCN and C2D 



show fractional abundance either boosted by up 
to four orders of magnitude as compared to the 
Models NS1-NS3 or stay unchanged, respectively. 
In the case of the high-z, the formation of DCN is 
dominated by the reaction D + HCN which is inef- 
ficient in our models for normal spirals because of 
the deficiency of atomic deuterium. In the models 
for normal spirals, the fractional abundance of D 
is indeed five orders of magnitude lower than in 
the models for high-z cores. 

3.2. D/H abundance ratios sensitivity 

As seen in Figs. 2] to [51 the D-fractionation 
calculated from the chemical abundances of all 
species do not show similar D/H ratio val- 
ues, regardless of the environment. In pro- 
tostellar galactic cores, such as IRAS 16293- 
2422, such behaviour is a l ready well-known (e.g . 



1995t ICeccarelli etTal] Il998t 



van Dishocck et al. 

Parise et al.l 12004 . Authors indeed observed 
NH 2 D/NH 3 =0.1 whereas they obtained DCO+ /HCO+=0.009. 
They interpret the higher D/H ratios they ob- 
tained as coming from a cooler and more extended 
gas located in the envelope around the source and 
not coming from the hot gas in the core. In our 
study, we can see (Fig. [4]) that the highest ratios 
are indeed obtained, on average, for the galactic 
cases (i.e. Models NS), for Model NSi which has 
got the smallest density (see Table [5]). 

In roughly all the models, there are three 
groups of molecules which converge to similar 
D/H values: the group including DCO+/HCO+, 
N 2 D+/N 2 H+ and HDCO/H2CO; the group con- 
taining NH2D/NH3 and D 2 CO/H 2 CO and the 
group containing CH 2 DOH/CH 3 OD, CH 2 DCN/CH 3 CN, 
CH3OD/CH3OH. D/H however varies by several 
orders of magnitude amongst those groups. With 
the current limited D/H dataset on galactic and 
extragalactic sources, confirmation of such results 
is not possible. 

The H2D + /H3 + does not vary as much as the 
other D/H ratios from a galaxy type to another. 
On average the largest discrepancy between D/H 
ratios group estimates is shown to be for the SB+ 
case where the spread between all the D /H values 
is the largest (see Figj5|). 

Despite these uncertainties and for providing 
future observers with D/H values, steady-state 
D/H ratios for various molecules are listed in Ta- 
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bleH 

3.3. Predictions for observers 

In Table [U we list the D-species that should be 
likely detectable in external galaxies as based on 
the values obtained for their chemical abundance. 
We define the limit of detectability as [n(X)/n ff ] 
= 1 x 10~ 12 , as is typical for dense gas molecules 
in the Milky Way. Of course, detectability will 
also depend on the excitation conditions as well 
as the geom etry of the source. A ccording to the 
estimates of Lintott et al however these 

abundances may be sufficiently large that unre- 
solved active galaxies should be likely detectable 
in these species even at high redshift. 

Of particular interest are HDO, DCN and 
HDCO which are very abundant ([n(X)/n#]> 
10~ 10 ) in all the investigated environments. These 
species are ideal molecules to be used for a first 
observational campaign aiming at validating the 
model predictions described here. 

A second category of interesting species in- 
cludes DCO+, DC 3 N, DNC and N 2 D+ as they 
seem to arise only from one gas component in a 
single galaxy type: SB 3 + (for DCO+), SBi (for 
DC 3 N), NS 3 (for DNC) and SB 3 + (for N 2 D+). 
In a similar vein, C 2 D and HDCS are predicted 
to be excessively reactive to the presence of cos- 
mic rays ionisation whereas they survive easily in 
FUV-enhanced environments. HDS is more sensi- 
tive to FUV changes than C 2 D and HDCS and it 
survives to both FUV and cosmic ray ionisation as 
long as the gas is dense. Here, HDS may thus be a 
good discriminant of the different gas components 
in SB and SB+ environments since it shows high 
abundances for the densest gas component. How- 
ever, as soon as the metallicity drops, it becomes 
undetectable whereas HDCS and C 2 D abundances 
are still high. 

The third category of species contains D 2 CO, 
CH 2 DOH and CH 3 OD. These molecules do not 
survive in high-redshift nor in low-metallicity 
galaxies. The fact that these species survive 
nonetheless in SB sources may reflect the poor 
influence of FUV photons on highly embedded 
material. 

Species such as CH 2 DCN and H 2 D + are likely 
undetectable, whatever the conditions since their 
chemical abundances are low. 



Finally, DCN, HDO, HDCO and to a lesser ex- 
tent D 2 CO, DCO+ and NH 2 D seem also to be 
good candidates i.e. abundant enough at various 
redshifts. 

Of course, our predictions are qualitative in 
nature since based only on an abundance crite- 
rion. For more quantitative predictions, radia- 
tive transfer models, as well as the knowledge of 
the source geometry, are needed. DCN, HDCO, 
D 2 CO, DCO+ and NH 2 D expected line inten- 
sities can not be derived due to the lack of colli- 
sional rates. For HDO however, we have been able 
to derive very rough estimates of line intensities 
using RADEX0 developed by Ivan der Tak et al 



(<2007l) and using the collisional rates from the Lei- 
den Atomic and Molecular Database (LAMDAE]). 
We have assumed a plane-parallel geometry, used 
Model SBi HDO fractional abundances for de- 
scribing the local Universe abundances, and Model 
High-zi HDO fractional abundances for the z> 1 
case. We have converted the predicted HDO frac- 
tional abundances from both models into column 
densities as described in Sect. 0](see the footnote), 
using the A v values listed in Table |U The HDO 
column densities obtained have then been used in 
RADEX, in addition to the gas density and the ki- 
netic temperature as listed in Table [2j A FWHM 
of 50 kms -1 has been assumed, as seen typi- 
cally in extragalactic molecu l ar lin e measurements 
(e.g. iMauersbereer et all 119951 llsrael fc Baasf 



20031: iMartm et all 120061: iBavet et all 120091 and 
ladro et al.l l2010h . For Models SBi, we have 
obtained values ranging from 1.8xl0 -4 Kkms -1 
(for the HDO(5 3j2 -6i j 5) line emitting at 356 GHz) 
up to 9.5xl0 3 Kkms" 1 (for the HDO(2 , 2 -li,i) 
transition emitting at 490 GHz). Beam and dis- 
tance dilution correction fact ors have to be a pplied 
to these values. Following iLintott et al ([2005) 
formula (see the Appendix), we have assumed a 
distance of 3 Mpc (such as M 82) as well as a 
beam size of 12 arcsec as listed for ALMA-band 
?H (i.e. at 490 GHz). With such assumptions, the 
HDO(2 , 2 -li,i) line intensity can be reproduced 
by ~ 10 5 — 10 8 cores, depending on the source 
size considered. If we perform the same calcu- 
lation but for the high redshift case (i.e. HDO 



1 See http://www.sron.rug.nl/~vdtak/radex/radex.php 



2 See http://www.strw.leidenuniv.nl/~moldata/ 
3 See 



primary beam size values in 
http: / /www.eso.org/sci/facilities / alma/observing/specifications / 
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fractional abundances from Model High-z-i), fo r 



a galaxy located at z=4.7 (jCarilli et al.l 12002f t . 



and increasing the ALMA beam to 56 arcsec 
(since the frequency line drops to ALMA-band 
1), we thus need ~ 10 9 cores to reproduce the 
HDO(lo,i-Oo,o) line emitting at 465 GHz, which 
shows the strongest RADEX line intensity esti- 
mate of 2.09X10 1 Kkms" 1 . In this case, HDO 
may not, of course, be detectable. 

4. Discussion and conclusions 

Although we do not aim at modelling any spe- 
cific source, we have qualitatively compared our 
models with extr agalactic D-species observa tions 



cept Model SBi+ can reproduce the observed up- 
pe r limit of N(N 2 D+) < 9.4 x 10 10 cm" 2 reported 



published so far. iMauersberger et al.l (|1995[ ) ob 



tained upper limits for the emission of the DCN 
J=2-l lin e toward NG C 253 and IC 342. More 
recently, iMartm et all (|2006l ) observed DCO+, 
DCN, DNC and N 2 D+ in NGC 253. When we 
convert these observations and our predicted frac- 
tional abundances into approximate column den- 
sitiet@, we find a general good agreement between 
the two sets within a factor of < 2-3, even if one 
single model is not able to reproduce consistently 
all the data. As already shown in iBavet et al 



(l2009h . several gas components are needed to 
account for the molecular emission observed in 
external galaxies. More precisely, in NGC 253 



(data from both IMauersberger et al.l 119951 and 



Martin et al1l2006l ). the Model SBi+ best repro- 
duces the observed DCO + column density (less 
than a factor of 3). This indicates that a cosmic 
ray-enhanced environment may be present in the 
nucleus of NGC 253. Model NS3 best reproduces 
the observed DCN column density in NGC 253, 
within a factor of 1. Finally, all models in our 
grid, except Model NSi, are able to reproduce 
the observation of DNC (N(DNC) < 7.2 x 10 10 
cm -2 ). It is doubtful whether we can compare 
the observed and the predicted column densities 
for N 2 D + because it is expected that doubly- 
deuterated species, not included in our current 
chemical network, might provide important for- 
mation and destructi on chemical routes for N 2 D + 
([Roberts et al.l 120031 ). With such restriction in 



mind, it appears however that all the models ex- 



4 We performed such conversion by simply multiplying our 
predicted fractional abundance by the column density of 
hydrogen at the relevant visual extinction. 



in IMartm et al.l (|2006h . For IC 342, only one D- 
species has been observed (i.e. DCN) and Model 
NS3 gives the closest predicted column densities 
to the observed values (factor of 11). This factor 
is higher than in the case of NGC 253 because 
the DCN detection suffers from larger uncert ain- 
ties in IC 342 (see IMauersberger et al.l Il995l ) as 
compared to the detection obtained in NGC 253. 

Due to the lack of additional observations in 
extragalactic environments, we are not able to 
further refine our predictions. Nevertheless, our 
simple approach offers some significant insights 
on the major trends of the deuterium chemistry 
when subject to various physical conditions. Our 
models can also potentially help discriminate cold 
and dense star-forming gas between SB, SB+, low- 
metallicity, etc activity. These predictions may be 
of particular interest when looking at extragalactic 
very dense high-mass star-forming regions where 
CO is expected to be depleted onto grains (simi- 
larly to what it is seen in our own Galaxy), and 
where thus D-species drive the chemistry. 
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Table 1: Standard model parameters (See Model NS3). 



P sx cLin6 1 6r 


Symbol 


Typical Milky 
Way values 


Collapse Mode a 


B 


0.1 


Initial number density (phase 1) 


rti 


300 H cirr 3 


Final number density (phase 1 and 2) 


n f 


lx 10 7 cm~ 3 


Temperature (phase 1) 


Ti 


10 K 


Temperature (phase 2) 


T 2 


300 K 


External UV radiation intensity 


I 


1 Habing 


Cosmic ray ionization rate 


c 


1.3xl0~ 17 s- 1 


Visual extinction 




580.6 mag 


Gas: dust ratio 


d 


100 


H2 formation rate coefficient 


R 


l.OxlO- 17 x VT cm 3 s -1 


Metallicity 


Z Q 


solar values, see Table [3] 



a : The collapse is treated as a 'modified free-fall', as defined bv lRawlings et all j|n)92|. and the parameter B 
is introduced to allow for collapse in free-fall B =1.0, or somewhat more slowly (B <1.0) if gas or magnetic 
pressure resists the collapse. 



Table 2: Initial elemental abundance ratios with 
respect to the total H-nuclei number, used as the 
standard values labelled ("ST") in Table M The 
initial elemental abundance ra tios correspond here 
to a compilation of v a lues f ro mlSembach &: Savage! 
(Il996h: ISofia et al.1 (fl997h : iMever et al.1 (Il998h: 
Snow et al. I (l2002h : iKnauth et all (120031 7 



ST 


C/H 


1.4x10" 


-4 


S/H 


1.4x10" 


-6 


O/H 


3.2x10" 


-4 


N/H 


6.5x10" 


-5 


He/H 


7.5x10" 


-2 


Mg/H 


5.1x10" 


-6 
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Table 3: Input parameters for the UCL_Chem models (see Sect. [5]). The abbreviation "ST" represents the 
standard values listed in Table Q]and Table [5J The abbreviation "CB-A, -B, -C, -D, -E" correspond to the 
values listed in Table [4] 



Model 


n f 


Size 


Metallicity 


Gas-to-dust 


Ini. Elem. 


c b 


Temp. 


I 


Comb. 




(cm 3 ) 


(PC) 


(z©) 


mass ratio 


Abund. ratios 


(Co) 


T 2 (K) 


(Ie) 


parameters 


Normal spiral 




















NSi 


10 5 


1.00 


1 


100 


ST 


1 


300 


1 


CB-A 


NS 2 


10 6 


0.15 


1 


100 


ST 


1 


300 


1 


CB-B 


NS 3 


10 7 


0.03 


1 


100 


ST 


1 


300 


1 


ST 


Starburst 




















SBi 


10 5 


1.00 


1 


100 


ST 


1 


500 


10 3 


CB-C 


SB 2 


10 6 


0.15 


1 


100 


ST 


1 


500 


10 3 


CB-D 


SB 3 


10 7 


0.03 


1 


100 


ST 


1 


500 


10 3 


CB-E 


Cosmic-Rays 




















Enhanced 




















SB! + 


10 5 


1.00 


1 


100 


ST 


100 


300 


1 


CB-A 


SB 2 + 


10 6 


0.15 


1 


100 


ST 


100 


300 


1 


CB-B 


SB 3 + 


10 7 


0.03 


1 


100 


ST 


100 


300 


1 


ST 


Low metallicity 




















Low- met i 


10 5 


1.00 


1/5 


500 


ST/5 


1 


500 


10 3 


CB-F 


Low-met2 


10 6 


0.15 


1/5 


500 


ST/5 


1 


500 


10 3 


CB-G 


Low-met3 


10 7 


0.03 


1/5 


500 


ST/5 


1 


500 


10 3 


CB-H 


High rcdshift 




















High-zi 


10 5 


1.00 


1/5 


500 


ST/5 


100 


500 


10 3 


CB-F 


High-z 2 


10 6 


0.15 


1/5 


500 


ST/5 


100 


500 


10 3 


CB-G 


High-z 3 


10 7 


0.03 


1/5 


500 


ST/5 


100 


500 


10 3 


CB-H 



a : zq — 1 corresponds to solar values of the elemental abundances ratios; b : expressed in units of (q = 
1.3xl0" 17 s" 1 . 
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Table 4: Combinations of parameters coupled with the metallicity and used in Table. [3] (see Sect. [5]). 





ST 


CB-A 


CB-B 


CB-C 


CB-D 


CB-E 


CB-F 


CB-G 


CB-H 


A v 


580.6 


194.9 


291.3 


194.9 


291.3 


580.6 


40.6 


59.9 


117.7 


(mag) 




















ratio of the number 




















densities of grains 


l.OxlCT 12 


l.OxlCT 12 


l.OxlO" 12 


l.OxlCT 12 


LOxlO" 12 


LOxlO" 12 


2.0xl0" 13 


2.0xl0" 13 


2.0xl0~ 13 


to hydrogen nuclei 




















H2 form." 


1.7xl0" ie 


1.7xlCT ie 


1.7xl0" ie 


2.2xl0" lb 


2.2xl0" lb 


2.2xl0" lb 


4.5xl0" 17 


4.5xl0" 17 


4.5xl0" 17 


rate coeff. 





















a : The H2 formation rate coefficient is computed at the beginning of the phase 2 of the UCL hot core model and it is expressed in cm 3 s 



Table 5: Deuterated detectable in extragalactic star forming regions and their corresponding D/H ratios (at a time=10 6 yrs) for five models 
representative of five types of galaxies (see Sect. [2J. The limit of detectability has been taken to be [n(X)/ns] = 1 x 10~ 12 , as is typical for 
dense gas in the Milky Way. Below this limit, we assumed the species as not detectable (symbol "-"). Otherwise, they are marked with the 



symbol "+". When the fractional abundance of a deuterated species is above [n(X)/n#] = 1 x 10 10 , the symbol "++" is used. 





Normal Spiral 






Starburst 










Cosmic- Rays 








Low 








High 
































Enhanced 








metallicity 








redshift 










NSi 


NS 2 


NS 3 


SBi 


SB 2 


SB 3 


SBi + 


SB 2 + 


SB3+ 


Low-mcti 


Low-mct 2 


Low-mct 3 


Higli-zi 


High-z 2 


High-z 3 


HDCO 


++ 


++ 


++ 




++ 


++ 




++ 




+ 




+ 


+ 


+ 




+ 


+ 




+ 


+ 




++ 




D 2 CO 


+ 


+ 


+ 




+ 


+ 




+ 
































DCN 


++ 


+ 


+ 




++ 


+ 




+ 




++ 




++ 


++ 


++ 




++ 


+ 




++ 


++ 




++ 




DNC 


+ 


+ 












































DC 3 N 










+ 






































DCO+ 




















+ 




























H 2 D+ 
















































HDO 


++ 


++ 


++ 




++ 


++ 




+ + 




++ 




++ 


++ 


++ 




++ 


++ 




++ 


++ 




++ 




J C 2 D 


++ 


++ 


++ 




++ 


++ 




+ + 












++ 




++ 


++ 




++ 


++ 




++ 




^HDCS 


++ 


++ 


++ 




++ 


++ 




+ + 












++ 




++ 


++ 




+ 


+ 




+ 




HDS 


+ 


+ 


+ 










+ 










+ 






















NH 2 D 


++ 


++ 


++ 




++ 


++ 




+ + 




+ 








+ 




+ 
















N 2 D+ 


























+ 






















CH 3 OD 


+ 


+ 


+ 




+ 


+ 




+ 
































CH 2 DOH 


+ 


+ 


+ 




+ 


+ 




+ 












+ 




















GH 2 DCN 
















































D 2 CO/H 2 CO 


5.1 x lCT" 


8.8 x 10"' 


6.2 x 10" 


6 


2.8 x 10 _o 


1.6 x 10" 


6 


4.0 x 10" 




1.2 x 10" 


4: 


4.2 x 10 -4 


2.7 x 10 _b 


3.1 x 10 


1; 


8.8 x 10 _ " 


3.7 x 10" 


ii 


5.5 x 10"' 


3.8 x 10" 


ii 


2.6 x 10" 


ll 


DCN/HCN 


1.4 x icr 4 


7.2 x 1(T 5 


2.0 x 10 


-4 


3.1 x 10~ 4 


7.6 x 10" 


-4 


4.2 x 10" 


-3 


1.6 x 10" 


-3 


8.7 x 10" 3 


2.0 x 10~ 3 


6.0 x 10" 


-4 


3.5 x 10" 3 


3.1 x 10" 


-3 


6.3 x 10" 4 


2.6 x 10- 


3 


8.2 x 10" 


3 


DNC/HNC 


4.8 x 1CT 5 


2.4 x 10 -5 


2.8 x 10" 




6.5 x 10" 6 


6.5 x 10" 


-6 


2.2 x 10" 




3.3 x 10" 


10 


2.3 x 10 _n 


5.4 x 10 -11 


1.9 x 10" 


-7 


5.9 x 10" 8 


8.8 x 10" 


-8 


9.6 x 10" 9 


3.2 x lo- 


9 


2.6 x 10" 


10 


NH2D/NH3 


2.7 x 10 -5 


1.3 x 10 -5 


1.4 x 10" 




3.3 x 10" 5 


2.4 x 10" 




2.6 x 10" 




5.0 x 10" 




1.5 x 10 -5 


3.0 x 10" 5 


1.3 x 10" 


-6 


1.5 x icr 6 


1.8 x 10" 


-6 


1.8 x 10 _0 


re x 10" 


6 


2.1 x 10" 


6 


DCO+/HCO+ 


3.9 x 10" 4 


2.0 x 10" 4 


2.1 x 10" 


-4 


3.0 x 10" 4 


2.4 x 10" 


-4 


4.2 x 10" 


-4 


6.5 x 10" 


-3 


2.2 x 10~ 3 


2.9 x 10" 4 


8.2 x 10" 




1.5 x 10" 4 


2.4 x 10" 


-4 


7.1 x 10" 4 


5.6 x 10" 


4 


5.4 x 10" 


4 


HDCO/H 2 CO 


1.0 x icr 3 


5.7 x 10" 4 


5.0 x 10 


-4 


,3.3 x 10" 3 


9.2 x 10" 


-4 


2.5 x 10" 


-4 


1.7 x 10" 


-2 


3.3 x 10" 2 


2.5 x 10 _1 


6.2 x 10" 


-3 


9.5 x 10~ 3 


1.4 x 10" 


-3 


7.5 x 10~ 3 


6.7 x 10- 


2 


1.4 x 10" 


1 


CH3OD/CH3OH 


1.1 x icr 3 


7.2 x 10" 4 


7.3 x 10" 


-4 


9.3 x 10" 4 


5.8 x 10" 


-4 


6.4 x 10" 


-4 


1.5 x 10" 


-2 


1.2 x 10~ 2 


3.6 x 10" 3 


8.0 x 10" 


-3 


6.7 x 10" 3 


9.1 x 10" 


-3 


4.6 x 10~ 3 


9.0 x 10" 


3 


8.3 x 10" 


2 


CH 2 DOH/CH 3 OH 


1.5 x 10 -3 


1.2 x 10 -3 


1.1 x 10" 


-3 


1.5 x 10 -3 


8.7 x 10" 


-4 


5.8 x 10" 


-4 


1.9 x 10" 


-2 


1.5 x 10 -2 


1.3 x 10 -2 


9.5 x 10" 


-2 


6.6 x 10" 2 


2.4 x 10" 


-2 


7.5 x icr 3 


3.0 x 10- 


2 


3.3 x 10" 


1 


CH 2 DCN/CH 3 CN 


7.2 x icr 4 


1.4 x 10" 3 


5.4 x 10" 


-4 


6.1 x 10~ 3 


1.0 x 10" 


- 3 


1.1 x 10" 


-4 


,3.1 x 10" 


-5 


2.9 x 10~ 5 


3.8 x 10" 5 


4.1 x 10" 


-3 


3.9 x 10" 3 


4.2 x 10" 


-4 


7.2 x 10" 3 


6.8 x 10- 


2 


5.3 x 10" 


2 


H 2 D+/H+ 


7.5 x 10-"' 


7.6 x lO" 5 


7.6 x 10" 




5.6 x 10" 5 


5.6 x 10" 




5.3 x 10" 


-6 


6.6 x 10" 




6.1 x 10~ 5 


5.6 x ID -6 


1.8 x 10" 




1.7 x 10" 5 


1.7 x 10" 




1.9 x 10~ 5 


1.1 x 10" 




4.7 x 10" 


6 


N 2 D+/N 2 H+ 


5.2 x 10~ 4 


2.1 x 10~ 4 


2.4 x 10" 


-4 


5.3 x 10" 4 


5.4 x 10" 


-4 


1.4 x 10" 


-3 


1.2 x 10" 


-3 


9.4 x 10~ 4 


1.9 x 10~ 4 


2.9 x 10" 


-4 
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Fig. 1. — Fractional abundances as a function of tinj^ of various D-species for the normal spiral case (top 
three plots) and the starburst case (bottom three plots). For each case, we have plotted the three gas 
components seen in Table[3]i.e the lowest density gas component (n(H2)=10 5 cm~ 3 , top), n(H2)=10 6 cm~ 3 
(middle) and the densest gas component (n(H2)=10 7 cm~ 3 , bottom). Fractional abundances are expressed 
with respect to the total number of hydrogen nuclei nu and the limit of detectability has been taken to be 
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2. — Fractional abundances as a function of time of various D-species for the cosmic-rays enhanced 
caption of Fig. Q}. 
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Fig. 3. — Fractional abundances as a function of tirnggOf various D-species for the low metallicity case 
three plots) and the high redshift case (bottom three plots) . See the caption of Fig. [TJ 
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Fig. 4. — D/H abundance ratios as a function of time for various D-soecies in the normal soiral case 
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Fig. 6. — D/H abundance ratios as a function of time for various D-species in the low metallicity case (top 



